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DECAY  OF  DROPS  BY  EVAPORATION 
IgriCce  I .  Kolodner 

1,   Introduction 

In  order  to  deter: line  the  radius  R(t)  of  an  evaporating 
drop  it  is  necessary  to  simultaneously  consider  the  concentra- 
tion of  vapor  c(r,t)  around  the  drop.   If  \<ie   assume  that  the 
teraperature  remains  constant  and  the  drop  remains  spherical, 
then  c(r,t)  and  R(t)  can  be  determined  by  solving  the  following 
problem: 

2 

(1)  c^^  +  -  c^  =  c^     for  r  >  majc  [0,R(t)]  ,   t  >  0 

c[K(t),t]  =  1  for  0  <  t  <  t^ 

ac^[R(t),t]  =  H(t)  for  0  <  t  <  t^ 

c(r,0)  =  0  for  r  >  1 

c(oo  ,t)  =  0  for  t  >  0 
R(0)  =  1 

c^(0,t)  =  0  for  t  >  t 

In  equations  (1)  all  quantities  are  in  the  dlmensionless 
form  given  in  [1]  (see  equations  (7) -(12)  of  [1]).   The 
constant  t   is  the  time  reqiiired  for  complete  evaporation,  and 
is  to  be  determined,  while  a  is  a  given  constant. 

The  main  results  of  [1]  are  that  over  most  of  the 
evaporation  time  the  evaporation  curve  is  approximately  given 

by 

(2)  R2(t)  -  1  -2at  -^  aJt      ' 

and  that  the  time  of  evaporation  t   is  approximately 
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(3)  t^^  =^[l._L_  ,7^^0(a)]   , 

3  'Stt 

Although  these  results  agree  well  with  experiment  and  with 
previous  approxiraato  results,  they  are  based  on  a  purely 
formal  approxiraation  method.   Therefore,  in  the  present  report 
the  problem  is  treated  rigorously  by  a  different  method. 
The  main  results  of  the  present  report  arc: 

(a)  The  problem  (1)  possesses  a  ujiique  solution  if 

■  0  <  a  <  1 ^   .38. 

1  +  JTW 

(b)  R(t)  is  the  limit  of  the  sequence  p  (t)  constructed 
in  section  3. 

(c)  For  every  n,  p^.^^lt)  <  R(t)  <  P2n+1^^^  '      '^^'■'^   ^^P^^^ 
and  loi,7cr  bounds  for  E(t)  arc  obtained. 


(d)   l-2a(t+2  i)  <  R2(t)  <  i.-i^t   . 
w  JL  —       —    a+1 


The  bounds  onR2(t)  in  (d)  arc  the  crudest  of  those  described 
in  (c).   Nevertheless,  the  lower  bound  coincides  with 
equation  (2),  the  result  of  [1],  and  if  the  lower  bound  on 
t^  in  (o)  is  expanded  for  small  a,  it  coincides  with 
equation  (3) . 

,  In  Figure  1  the  two  bounds  in  (d)  are  plotted  for 
c  =    .05,  and  an  improved  upper  bound  is  also  shovm.   A  large 
value  of  a  has  piirposcly  been  chosen  in  order  to  show  the 
differences  betx,reen  the  bounds.   In  a  practical  problem  in 
which  a  is  of  the  order  lo"^  the  various  bounds  could  not  bo 
distinguished  from  each  other,  and  therefore  either  of  them 
would  represent  the  solution  adequately. 
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The    dlffcroncG   between   the   two   bounds    in    (d)    is   less 

thssi  — -y(l  +  Lj.l — )    during   the    full   time    of  evaporation,   vjhero 
a      is   the   lov;er  bound   in   (c).      Similarly,    the   bovinds    in    (e) 


a 


-1 


a 


determine  t  with  a  relative  error  of  less  than  \  —^^r"  " 
o  L  2a 

For  a  water  droplet  evaporating  in  air  at  73  Pj  g-^-  2  ?4  10  -^ 
and  these  errors  are  .0071  and  ,0072  respectively. 

The  restriction  on  c  in  (a)  is  a  consequence  of  the 

method  of  solution.   In  fact,  for  0  <  a  < s^i  .57  the 

method  still  yields  a  solution  but  not  for  all  values  of  t  up 

to  t  ,   The  solution  so  obtained  would  have  to  be  extended  up 
o 

to  t   by  another  method.   Most  likely  the  solution  does 
exist  for -all  a  >  0,  but  the  present  method  vjill  not  apply 
for  all  a,  since  it  can  be  shown  to  fail  if  .90  <  a  <  I.15. 
In  the  corresponding  one-dimensional  problem  the  present 
method  does  yield  a  solution  for  all  a  >  0. 
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.  If  a  <  0,  the  vapor  condenses  and  the  drop  gro\<is,    so 
R  >  0.   Therefore,  Theorem  3  of  section  2  is  not  applicable, 
and  the  present  raothod  must  be  modified.   A  procedure  for  doing 
this  has  been  found,  and  will  be  reported  in  the  near  future. 
In  the  one-dimensional  case  the  exact  solution  is  l-aioim   for 
all  a  >  -1  and  t  lorefore  the  same  restriction  on  a  must  bo 
expected  in  the  spherical  case. 

In  section  2  we  derive  an  integral  representation  for 
tne  .  condon.:ation  depending  on  the  evaporation  curve  R(t), 
and  show  that  this  expression  satisfies  all  of  the  conclitlons 
Imposed  on  c(r,t),  equations  (1),  provided  that  R(t)  is  the 
proper  evaporation  curve.   T^(t)  is  determined  as  a  ,-olution 
of  a  certain  integro-dlf fcrential  equation.   In  section  3, 
the  equation   for  R(t)  is  solved  by  constructing  a  sequence 
of  fujictions  which  converge   to  the  solution.   In  section  [}. 
x/e  derive  simple  loirer  and  upper  bounds  for  R(t).   (see  above) 
Finally,  section  5  contains  some  cor.iments  concerning  tne 
method  and  its  applicability  to  other  problems  of  similar 
structure . 

The  proof  of  statc.ients  appearing  in  sections  2-1;  is  in 
some  cases  long  and  intricate.   We  therefore  eliminated  all 
those  proofs  from  the  t.xt,  in  order  not  to  obstruct  its 
clarity.   Thoy  will  be  published  in  a  future  issue  of 
Cormiumications  on  Pure  and  Applied  Mathematics,  IJew  York 
University.   In  the  same  spirit,  statements  of  theorems  have 
been  simplified  in  some  cases  so  as  to  olim.inate  certain 
"minor"  points. 

2 •   Integral  representation  for  c(r.t) 

In  place  of  c(r,t)  it  is  more  convenient  to  consider 
the  function. 

^-1-)  u(r,t)  =   rc(r,t)   . 

I^cm  (1)  it  follows  that  u  satisfies  the  following  conditions: 
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(>)         Upjo   "^  ^t      '  ^   ^  ^^^^    (0,R(t))     ,       t   >    0 

(6)         u(R(t),t)    =.  R{t)       ,  0   <    t    <    t 


o 


(7)         au   (R(t),t)    =  R(t)R{t)  +0.   ,  0  <   t  < 


r 


u. 


O 


(8)  u{v,0)    =    ')      ,  r  >  1    ,      u(0D,t)    =   0 

(9)  Ti{0)    =   1      , 

(10)  lim  H  exists   and   is   bounded   for   t  >  t 
r— >0^'  ° 


ru     -  u 
(11)         lim  — ^ —  =   0    .      t  >  t 
r— ^0 


V —  =   0   *       - 
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vie  first  consider  the  problem  in  the  range  0  <  t  <  t  and 
consequently  disregard  the  conditions  (10)  and  (11)  x-zhlch 
a-PPly  only  for  t  >  t  , 

Let  U(r,t,^,<3-)  be  the  fundamental  solution  of  the  heat 
equation 

(12)  y  =  —^rrtr  exp  (-(-^-^J^)        . 

2/it(t-ar)    V  2.[t^  / 

If  the  solution  to  cur  problem  were  known  to  exist  and  if  the 
boundary  f-onction  R(t)  v/ore  known,  then  u(r,t)  could  be 
obtained  from  the  well-lcnown  representation  of  solutions  of 
the  heat  equation  by  line  integrals  involving  boundary  values, 
viz; 

(13)  u(r,t)  =  j>     u(^,(r)Ud^+rij||-u||]d(r  , 

see  e.g.,  [2],  vol.  HI,  p.  31.I,  formula  (35).   Substituting 
the  data  (6),  (?)  and  (0)  into  (I3)  wo  obtain  the  formal 
expression 
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v/licro 

(15)  z  =  z(P,cr)  =  --JiZ}      , 

zft-ar 

(16)  Y7(z)  =  -i  exp  (-z^)   , 

and  r  >  R ( t ) ,  t  >  0 . 

It  is  clear  that  c;:prcssion  (llf)  with  arbitrary  n(t) 
cannot  satisfy  all  the  conditions  ImposGd  on  u,  since  then 
the  boundary  would  be  proscribed,  and  by  a  vroll-knovm  theorem 
a  solution  of  the  heat  equation  is  uniquely  determined  by 
specifying  one  boujidary  condition  only,  and  not  two  as  in  the 
present  case,  sec  equations  (6)  and  (?)•   One  uaj   to  proceed 
would  be  to  determine  R,(t)  by  requiring  that  either  one  of 
(6)  or  (7)  be  satisfied.   This  apparently  simple  procedure 
presents  a  number  of  difricultios ,   Firstly,  it  is  not  clear 
whether  equation  (li|)  or  its  derivatives  coixld  be  used  for 
r  =  H(t),  since  its  value  on  the  boundary  (if  it  exists  at 
all)  may  not  be  equal  to  the  limit  from  the  Interior, 
Secondly,  if  the  first  difficulty  could  be  overcome,  the 
resulting  equation  for  R(t)  would  not  be  easy  to  solve. 
Lastly,  even  if  R(t)  could  be  so  determined,  it  would  again 
be  difficult  to  show  tha,t  the  second  of  the  two  conditions  (6) 
or  (7)  is  then  satisfied. 

The  object  of  the  remainder  of  this  section  is  to 
determ.ine  an  equation  for  R(t)  ^.w!   to  show  that  vjith  the 
proper  R(t)   (l.'-|-)  actually  represents  the  solution  to  our 
problem,  in  the  evaporation  stage.   To  this  effect,  wo  consider 
an  arbitrary  smooth  fimctlon  p(t)  which  decreases  monotonically 
from  unity, 

(17)  p(0)  =  1   ,    p(t)  <  0 
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(p  need  not  be  finite  cvoi>yT-;hGro )  .   We  form  the  fimction 
u'°(r,t)  by  replacing  R(<cr)  by  p((r)  in  (lii)  and  (l5).   That  is, 

t 

(18)  up(r,t)  =1  rr^^^ 


{i+c~^)n{cr)  0 
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iii£lolqiL±i)Y.(z)d(r 

t  -  <r      y 


This  ftinction  is  defined  here  for  both  r  >  p(t)  and  r  <  p(t), 
and  for  t  >  0.   Observing  that 


1  P(^)  y      1+  o(CT')o(ar) 


Y|(z)  =  r/|(z)^+  («/t"-Crr^(z))^   , 


0  <  cr<  t    , 

wc  obtain 

z(r,t) 
(19)   uP(r,t}  =  r  r      y^(r)dr-jT  Vi(z(r,0)) 

^z(r,0) 


1     C     P(O')  P((r)     ;   X  .^ 


Upon  differentiating  (19)  we  find 

z(r,t) 

t 


(20)   uP(r,t)  =  r   .    Kj(t;)d2r  -  -^r^(z(r,0)) 


■z(r,0) 


+ 


-L  f  p(g>^)  p(g')  ^  n(-)d<r 

2a  J     t^^O^^   -n^-)^^*J 


Since 
(21) 


z(r,t) 


+00   if  r  >  p(t) 


-CD   if  r  <  p(  t) 


uP,  u^  are  certainly  discontinuous  across  the  line  r  =  p(t). 
We  denote  the  limits  of  u^,  u^   as  r  — 2>-p(t)  from  the  right 

(loft)  byuP(p^,t),  uP(p'',t)  (uP(p",t),  uP(p~,t)). 
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Theoreml: 

a)  u'  is  a  aolution  Ox"  the  hear  equation  for  r  <  p(t) 
and  r  >  p( t),  and  t  >  0. 

b)  Expressions  (I9)  and  (20)  have  a  meaning  for  r  =  p(t.), 
which  we  denote  by  uP(p,t),  <(p,t). 

c)  uP(p-,t)    =  uP(p,t)    +  ip(t) 


r 


i^(p".t)   =  uP(p,t)    i  J(a"-p(t)p(t)+  1). 
d)      u^  -*0  as   t  -»  0    (r  7^  0),    axid  also  as    irl 


CO 


Ho  observe  that  part  c)  ccucerning  u^  follows  easily  if  wc 
assume. the  continuity  of  the  last  term  in  equation  ( I9 ) . 
Indeed, 


±co 
+ 


uP(p-  ,t)  =  p  r      »^('t)dr  +  A 
z{p,(.) 

0 

u'^(p,t)  =  p     r      .       Y^(2r)d^+  A  ,  since  z(r,t)  =  0, 
zip>0)  for  r  =  p(t) 


and 


uP(p-,t)  -uP(p,t)  =p  r    ^(tjdr--.  t  ip(t)  . 

^0 


A  stands  for  the  last  two  terms  in  (I9)  evaluated  at  r  =  o(t) 
Vie   now  determine  the  solution  of  the  equation  (22) 


(22)  uPfn" 


u^(p",t)  =  0   ,  p(o)  =  1 

Vihlch    solution  wo    deno-L-e   by  R(t). 

Jlicorem_2:      Equation    (22)    has   a  unique    solution  R(t)    for 
'■^1*1  ^o>    and   this    solution  is   smooth.      Furthermore, 
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R<0,   0<t<t   ,   R(t)=0   . 


This  solution  is  constructed  in  section  3» 
Theorem  3: 

u^(p,t)  =  0  for   r  <  R(t)  ,   t  >  0   . 

Wc  now  can  prove  the  main 

Theorem:   u(r,t)  =  u  (r,t)  for  r  >  R(t)  satisfies  conditions 
(5)  to  (9). 

Proof;   By -Theorem  2,  R  <  0,  hence  Theorem  1  applies  to  vP 
TTlth  p  =  R,  and.  therefore  (5)  and  (8)  arc  satisfied.   By. 
Theorem  3,  u"  (r,t)  =0  for  r  <  R(t),  hence  in  particular, 
u^(R",t)  =  0,  \'Jo   now  have  from  Tlioorem  1, 

u^(R'^,t)  =  u^(R",t)  +R(t)  =  R(t) 

u^(R'',t)  =  u5(R",t)  +  (a"^v(t)F.(t) +  1)  =  a'V.(t)R(t)  +1  , 

hence  conditions  (6)  and  (7)  are  satisfied.   This  completes 
the  proof.   Obviously,  'S^{r,t)    is  identical  with  {lb.)    for 
r  >  R(t) . 

To  construct  the  solution  for  t  >  t   as  vjell,  wc  consider 
the  function 

-00 


(23)   v(r,t)  =r  r   .   >^(r)dr-  ft  rj(  z(r,0)  ) 

'-'^/■v,  .■:^ 


z(r,0) 


rain  (t,t^) 


where    -oo    <  r  <   oo  ,    t  >   0.      V/o    show  that 

(2I4.)      u  =  r+v(r,t)  -v(-r,t)    ,      r  >  max   (0,R(t)) 
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satisfies  all  the  conditions  (5)  tlirough  (11).   It  follows 

from  Theorem  1  that  v  is  a  continuous  function  of  r  and  t  for 

0  <  t  <  t  ,   It  can  be  shov;n  that  v  is  an  analytic  fvmction  of 
—  o 

r  for  t  >  t   and  a  continuous  function  of  r  and  t  for  t  >  t  , 
o  —  o 

see  [2],  ch.  XXIX.   Clearly,  for  t  <  t^, 


.R 


v(r,t)  =  Tx^'(r,t)  =0   ,      for  r  <  R(t) 


R 


v(r,t)  =  u(r,t)  -r  =  u-r   ,  for  r  >  R(t) 


Hence 


v(-r,t)  =  0   , 


for  r  >  -:;(t)   , 


and  therefore  for  r  >  U{t)    since,  R(t)  >  0.   ConseqUGntly, 


u  =  u^Mr,t)   , 


for  r  >  R(t)  ,  ■fc  <  "^Q  5 


as  required.   For  t  >  t  we  use  the  analyticity  property  of  u. 

Since  u  is  an  odd  function  of  r,  —  is  even  and  analytic. 

ru^  -  u 

Hence  —  exists  as  r  — ^  0  and  i:;^)      -   — ^ —  =  0  as  r  — s>  0. 

"       r 
This  completes  the  proof. 

Explicitly,  the  final  formula  for  c(r,t)  is,  after  some 

transformations, 

-00 


(25)   c(r,t)  =  1  + 


J    •    (i  -  ^z^yi{v)dv 


-   z(r,0) 


min  (t,t^) 


2ar 

^0 
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R(cr)  R(cr) 


r^(z)d(r 


min  (tjt^) 
L^z(-r,0)  0 


-oo 


Here  t  >  0,  r  >  max  (0,R(t)),  and  R(t)  is  determined  by 
equation  (22).   The  second  bracket  Is  identically  zero  for 


t  1  ^o- 
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3 ,   Determination  of  the  evaporation  curve 

The   evaporation   ciirvc    is   determined  by.  equation   (22)    of 
the   preceding   section.      In  view  of   Ttieorem  1,    this   oqtiation  is 

(26)  uP.p,t)    =  |"a"^-p(t)p(t)  +1;.;       ,      p(0)    =   1 

or  explicitly 

■  0 

(27)  ^fa'^p(t)p(t)  +  ij  =  f  y|{'^)d'if-  -i-y|(z(p,0)) 

0_  /  -   r^\  2.PC 


''z(p,0)  2^ 


0 

After  some  transformations  this  can  be  reduced  to  a  simple] 
form 


(28) 


where 


rpp  =  G{p) 


2(P,0) 

(29)   G(p)  =  -2a  r     (1--^)  yi{t)o.t 

J-co       ^ 


0 

This  equation  is  best  handled  by  the  method  of  iterations 
V/e  first  define  an  auxiliary  function  x^(t,a,k)  by 

(30)  w^ft,a,k)    =  k^  -2c[(t-a)  +^{{1  -JK   )] 

where 


2~^       7  \2 


(31)      a>0,    0<k<l,    0<t<  t(a,k)    -  f  -  -i  +    rl^L2l+l 


+   (  i!  ■  ;?  i 


T.i ' 


'  \ 
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w(t,a,k)  passes  through  '^he  point  with  coordinates  (a,k) 
and  vanishes  for  t  =  t(a,k),   Furtherraore , 

(32)        ww'  =  -a(l+  -^)  =  G(l)  ,   t  >  0   . 


V/e  now  form  the  sequence  of  functions  p^^lt)  and  the  sequence 


of  numbers  a  ,  k  as  follows: 
n   n 


{53)  Po'vt)  =  w(x,0,l)  ,   0  <  t  <  t(0,l)  =  a^ 


t 


(33")    P2n+l^*^  =  ]^^'''  J     G(p2^)dt)  ,   C  1  t  <  a.^^^^ 


'0 

^2n+l  "  P2n+1^^2n+l^  ' 


v;h 


lere  a^   .,  is  the  first  zero  of  p„   n(t), 
2n+l  '^2n+l'   ' 


t 


(33"")  rP2nf2(^>  =2(1+  J'  G(p2^^^)dt  ,   0  <  t  <  a^^.^^ 


-  w(  t ,  a2^^4.;L »  ^2n+2  ^ '  ^2n+l~*-*  ^  ^2n+l '  ^'2n+2  ^  "^2n+2 


where  k2n+2  =  P2n+2^^^2n+l^ 

1 


Tlieorem:   If  a  < 


1  +J1j7% 


a)  The  seouence  a   converges  to  a  finite  mx'.iber  t^ 
'        -       n       '^  o 

b)  The  sequence  k   converges  to  zero, 

c)  For  any  t  <  t   there  is  an  N  such  that  for  n  >  N 
each  p^„  is  defined  for  all  0  <  t  <  t  and  the  sequence 
p   converges  to  a  fuiiction  R(t)  which  is  continuously 
differentiable  (except  at  t  =  0).  R(t)  is  the 
unique  solution  of  (23).   Furthermore, 

ILn  R(t)  :--^  0  ,   R(t)  <  0   . 
t— ^t„ 


(•;. 


I  \  -■ 
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d)  For  any  n, 

PPn  ^  P2n+2  ^  ^  ^  P2n+3  -  P2n+1 

P2n  "^  P2n+2  <  ^-  <  P2n+3  <  P2n+1 

whenever  any  o£   tv;o  compared  fionctions  are 
simultaneously  defined. 

In  viovj  o?  this  theorem,  the  evaporation  curve  is  constructed 
ever  the  full  period  of  evaporation  by  means  of  a  single 
sequence  of  functions.   Furthermore,  in  viev;  of  part  d),  any 
two  consecutive  iterates  yield  lower  and  upper  bounds  for  the 
solution  and  its  derivative, 

k-      E'stimatcs  for  R(t) 

Using  part  d)  of  the  theorem  of  the  preceding  section, 
V7e  have,  for  any  n, 

^^^^         PPn  ^  ^  -=  P2n-M  '      P2n  ^  ^  ^  P2n+1   ' 

Here,  the  first  part  of  the  two  inequalities  holds  for 
0  <  t  <  a^  ,  -whereas  the  second  part  holds  for  C  <  t  <  a-^ 
In  particular,  if  a  is  snail,  one  gets  simple  estimates  which 
are  sufficiently  sharp  for  all  practical  purposes  by  taking 
n  =  0  in  (3I4-).   Thus  we  get,  with  w  =  w(t,0,l),  (see 
equation  ( 30) ) 

(35)  R  >  w  ,   R  >  w  ,   0  <  t  <  a 


G(w) 


l2(l+  r    G(w)dt) 
J      ^0 


0  <  t  <  a^  <  a 
—   —  1     o 


^'      I 


III-. 

The  nurabors  Up  ,  and  in  particular  a  yield  lower  boionds 
for  t  .   None  of  the  cabcve  estimates  yields  an  upper  bound 
for  t  .   Such  a  bound,  together  v/ith  sharp  upper  bound  for  R(t) 
(though  not  for  R(t))  is  obtained  in  a  strikingly  simple  way 
directly  from  the  formulation  of  the  problem.   As  is .well  knovm, 
sec  e.g.  [2],  ch.  XXIX,  for  any  regular  solution  u(r,t)  of  the 
heat  equation  in  a  closed  domain  D, 


(36)  y (u^dt  +  udr)  =  0    . 

Wo  apply  this  identity  to  a  strip  bounded  by  the  derivatives 
t  =  t',  t  =  Oj  and  the  evaporation  curve,  obtaining 

■00         ■  t  '  .  CO 

(27)   j     u(r,t')dr=-r  (u^(R, t )+u(R, t)R( t ) )dt  +  T  u(r , 0)dr 
R(t')  ^  ^1 


_   -!- 1   a+1  ,  _,2  /,'v   TV    ^.  ^.i  ^. 
" 2a  (^(t)-l),   '-<^-  ±  ^ 


Q 


Now,  it  can  bo  shown  by  using  the  minimLim-maxlmum  principle  for 
solutions  of  the  heat  equation,  see  [3]?  that  n  >  0  fo.r  all  t 
and  r   .   Hence  wo  get 

(38)  r2  <  1-^t  ,   0  <  t  <  t,^   , 

and  since  R  >  0, 

These  bounds  have  been  discussed  in  the  introduction. 


This  is  obvious  on  physical  grounds,  siiice  n  represents  a 
density  concentration. 


■11' 


r-..r>p 


\'.i. 


O.-^n    \. 
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5,  Remarks  on   the  method  of  solution 

The  main  advantage  of  the  method  used  here  to  solve 
the  equation  for  R(t)  is  that  it  yields  in  a  single  sequence 
a  solution  over  the  full  range  of  the  evaporation  time,  and • 
thus  permits  an  estimate  of  the  evaporation  tine,   jfloroover, 
this  sequence  is  alternating  so  that  any  two  successive 
Iterates  yield  upper  and  lower  bounds  for  the  solution. 

The  final  expression  for  the  concentration,  equation 
(25)  could  bo  derived  by  evaluating  first  the  Laplace  trans- 
form of  the  concentration  and  inverting.   Since  to  determine 
the  Laplace  transform  one  has  to  solve  a  second  oj.^dor 
ordinary  differential  equation  with  constant  coeff ^".cicnts, 
and  since  the  inversion  can  be  carried  out  with  the  help  of 
existing  tables,  such  a  procedure  would  seem  at  first  sinplor. 
and,  espcclalljr  more  direct  than  the  cnc  used  here.   Actually 
one  would  thus  obtain  only  a  formal  expression  which  it  is 
difficult  to  handle.   We  would  have  at  our  disposal  the 
formula  (25)  from  which  K(t)  would  have  to  bo  determined  in 
some  way.   The  matter  v/ould  be  complicated  by  the  appearance 
of  the  second  bracket  in  (25)  which  as  we  knovj  novi  is 
identically  zer'o  during  the  evaporation  time,  but  which  x^^ould 
persist  in  an  equation  determining  R(t)  if  we  use  equation 
(25)  as  the  starting  point.   This  may  be  vjhy  no  adequate 
treatment  of  this  apparently  simple  and  fundamental  problem- 
has  appeared  until  no^^T,   The  method  used  hero  avoids  the 
difficulty  and  leads  to  equation  (25)  as  an  end-product  and 
not  as  a  starting  point. 

We  note,  without  fuj?ther  elucidation,  that  the  method 
illustrated  here  can  and  has  been  successfully  used  in 
handling  various  problems  of  similar  structure,  such  as 
problems  of  change  of  phase  with  unknown  (free)  bounidarlcs . 
Sec  [UJ. 


w- 
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